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ABSTRACT 

We present a new way of looking at the very long term evolution of GRBs in which the disk of material 
surrounding the putative black hole powering the GRB jet modulates the mass flow, and hence the efficacy of 
the process that extracts rotational energy from the black hole and inner accretion disk. The pre-Swift paradigm 
of achromatic, shallow-to-steep “breaks” in the long term GRB light curves has not been borne out by detailed 
Swift data amassed in the past several years. We argue that, given the initial existence of a fall-back disk 
near the progenitor, an unavoidable consequence will be the formation of an “external disk” whose outer edge 
continually moves to larger radii due to angular momentum transport and lack of a confining torque. The mass 
reservoir at large radii moves outward with time and gives a natural power law decay to the GRB light curves. 

In this model, the different canonical power law decay segments in the GRB identified by Zhang et al. and 
Nousek et al. represent different physical states of the accretion disk. We identify a physical disk state with 
each power law segment. 

Subject headings: gamma rays: bursts 


1. INTRODUCTION 

Gamma-Ray Bursts (GRBs) are thought to be energetic 
events heralding the formation of black holes (BHs): the ex- 
plosion of a massive star at cosmological distances gives rise 
to a long GRB, or two neutron stars coalesce to make a short 
GRB. (For a recent review of GRBs, see Gehrels, Ramirez- 
Ruiz, & Fox 2009.) To see a GRB at all, our line-of-sight 
must lie close to the BH angular momentum vector. A jet 
is launched toward the observer as the BH forms. The ini- 
tial Lorentz factor of material in the jet ~ 10 2 - 10 3 is large 
enough so that the observer only sees the central portion of the 
jet, due to relativistic beaming. The transition from “spher- 
ical” to “jet-like” expansion of the GRB ejecta is the basis 
of the current framework within which the analysis of GRB 
light curves is carried out (Rhoads 1997, 1999; Sari, Piran, & 
Halpern 1999). Parts of the jet further away from the sym- 
metry axis are initially invisible to the observer, hence the ex- 
pansion is indistinguishable from that of of a relativistically 
expanding sphere. Eventually interaction with the circumstel- 
lar medium reduces the Lorentz factor until the edges of the 
jet are visible, after which the decay of the GRB afterglow 
light curves becomes more rapid. This fundamental property 
of the decay relies only on special relativity and is not spe- 
cific to any one wavelength so these theoretical decays should 
have a similar, strong “break” from shallow to steep across all 
wavebands, hence the designation “achromatic break”. 

Frail et al. (2001) studied the decay properties of 17 pre- 
Swift GRBs with redshifts. By analyzing the multiwavelength 
decay properties, they were able to determine putative break 
times, which, when combined with distances and isotropic- 
equivalent GRB energies, allowed a correction to the true, 
beaming-corrected energy. Frail et al. (2001, see their Fig. 2) 
show' how a spread in isotropic energies between ~ 10 ,)2 and 
~ 10" 4 erg for their sample is considerably reduced down to 
a standard, geometry-corrected energy centered at ~ 3 x 1CP 0 
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erg. A detailed examination of the data entering into their 
study (see their Table 1 ) shows the sparse coverage of many 
of the light curves and underscores the difficulty in reliably as- 
signing break times 3 . An updated study along the same lines 
as Frail et al. (2001) by Bloom, Frail, & Kulkarni (2003) con- 
tains a somewhat larger sample — 28 GRBs with redshifts. 
We should emphasize that we do not argue against the need 
for beaming; clearly energies as large as ~ 10° 4 erg are prob- 
lematic, given that M e c 2 = 1.8 x 10 54 erg. It is important 
to distinguish, however, between a time-variable relativistic 
beaming, and geometric beaming. In this work we adopt a 
time-constant geometric beaming factor /beam = 3 x 10“ 3 to 
convert from isotropic energy to true energy. The question of 
interest to us is whether the very long term decay in the flux 
from GRBs, particular the X-ray flux, is due to the decelera- 
tion of a jet, or a decreasing mass supply rate onto the central 
engine. 

It was anticipated that Swift (Gehrels et al. 2004) would 
provide good examples of shallow-to-steep achromatic de- 
cays. After > 400 Swift GRBs, and now a total of > 170 
redshifts for all GRBs (41 for pre-Swift GRBs), this has 
not turned out to be the case. In spite of many GRBs 
with exquisite wavelength coverage, observed over many 
e — foldings of decay in flux, there is only one GRB in the 
Swift era perhaps showing an achromatic break, seen in X-ray 
and optical (GRB060206 - Curran et al. 2007) and a small 
number with candidate jet-breaks (Blustin et al. 2006, Stanek 
et al. 2007. Dai et al. 2007, Willingale et al. 2007, Kocevski 
& Butler 2008); the overwhelming majority do not show one 
(e.g.. Oates et al. 2007; Racusin et al. 2008; Liang et al. 
2008). Liang et al. (2008) present a thorough investigation of 
Swift XRT and optical jet break candidates. They note that "It 
is fair to conclude that we still have not found a textbook ver- 
sion of a jet break after many years of intense observational 
campaigns.” 

Zhang et al. (2006, see their Fig. 1; also Nousek et al. 
2006, see their Fig. 3) present a simple schematic for the de- 

3 The data used by Frail et al., as well as that for all GRBs, is viewable at 
http : / / grbiog . org/grblog . php 
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Fig. 1 . — A schematic representation of the X-ray decay of a Swift GRB 
as seen with XRT, based on the Zhang et al. (2006) idealization and GRB 
060729 from Willingale et al. (2007, see their Fig. 10). To account for cos- 
mological time dilation for a given GRB, the x — axis needs to be shortened 
by a factor 1 + z. 

caying GRB light curve as seen by the XRT on Swift. The 
decay is traditionally shown in logF 1 - logf. There are four 
basic power-law decay (F oc t~ n ) regimes: (I) o-i ~ 3 out 
to 10 2 — 10 3 s, (II) an — 0.5 out to 10 3 - 10 4 s, (III) 
am ~ 1.2 out to 10 4 — 10 5 s, and (IV) aiv — 2 at late 
times. The general picture given in Zhang et al. based on 
a handful of case studies from early Swift XRT results has 
been borne out by studies relying on a much larger sample 
(e.g., O’Brien et al. 2006, Ghisellini et al. 2009). In addition 
one also frequently sees flares within segment II (Burrows, et 
al. 2005; for possible explanations see Perna, Armitage, & 
Zhang 2006, Lazzati, Perna, & Begelman 2008), which are 
beyond the scope of this work. Figure 1 shows a schematic 
reproduction of the idealized Zhang et al. light curve, with 
realistic numbers for the XRT flux. The conversion from the 
XRT flux values to a beaming-corrected luminosity is given 
by Lxrt — 47ni| /beam i'xR.T* where dL is the luminosity 
distance. 

In this work we present a framework based on the long term 
evolution of a transient accretion disk formed out of some 
small fraction of the progenitor to account for segments II, 
III, and IV. (Segment I is too steep to be accounted for in our 
formalism.) The flatness of segment II arises from the tran- 
sient adjustment of a small amount of gas into an accretion 
disk, the segment III decay is consistent with the self-similar 
decay of a disk without a confining external torque, and the 
increase in slope in segment IV would be caused by the onset 
of significant disk evaporation (or by the late-time accretion 
from a super-Eddington slim disk). The consistency of the 
theoretical decay with segment III argues for a direct associa- 
tion, and the fact that it takes some time for an arbitrary initial 
mass profile to adjust into a disk leads to the transient interval 
(segment II). 

This work applies mainly to the commonly observed “long” 
GRBs which are believed to be due to the core collapse of a 
massive, high angular momentum progenitor, leading to the 
formation of black hole — the “collapsar” model (Woosley 
1993;MacFadyen& Woosley 1999; Woosley & Bloom 2006; 


Nagataki et al. 2007). There follows a period of electromag- 
netic extraction of energy from the BH which powers a colli- 
mated jet. For ~ 10 2 — 10 3 s following the phase of prompt 
GRB emission, ~ 0.1M© of bound matter from the disrupted 
progenitor star returns to the BH vicinity due to conservation 
of angular momentum. Kumar, Narayan, & Johnson (2008ab) 
attempt to explain the long-term Swift X-ray light curves in 
terms of various phases of an accretion model. They argue 
that the rapid accretion of the outer half of the core of the pro- 
genitor, possessing a density profile that drops steeply with ra- 
dius, can lead to the initial, rapid decay. For the flat segment 
(II) they discuss briefly the idea of having a small accretion 
disk viscosity, but ultimately discard this explanation in favor 
of prolonged fail-back of gas from the envelope of the progen- 
itor at larger radii, which they argue can supply matter to the 
central engine at a near-constant rate, or a rate that decreases 
slightly with time. Kumar et al. (2008a) also present analytic 
time dependent solutions for an accretion disk in an ADAF- 
like state (ADAF = advection dominated accretion flow) and 
show that, for a radial mass profile M(r) oc r s , the decay 
law would be am = 4(.s + l)/3. For s = 0, i.e., no wind 
mass-loss from the disk, the decay is close to the Zhang et 
al. (2006) value. Under the assumptions of core fall-back for 
segment I and envelope fail-back for segment II, Kumar et al. 
(2008b) place constraints on the progenitor core and envelope 
radii, their density profiles, and the pre-SN rotation rates of 
the core and envelope, normalized to their Keplerian values. 

In Section 2 we discuss the physics of accretion disks, both 
steady state and time dependent. Section 3 presents a scenario 
for the overall, long term evolution of the external accretion 
disk, section 4 presents a discussion, section 5 discusses the 
implications of the theory, and section 6 sums up. 

2. ACCRETION DISK PHYSICS 

Most of the previous discussion of GRB fall-back disks has 
failed to recognize one unavoidable aspect of their existence: 
namely, their outer edge moves outward with time, carrying 
with it mass and angular momentum. (Notable exceptions are 
Kumar et al. [2008ab].) One can make simple arguments 
about the perturbation accretion time scale near the central 
engine based on physical conditions there, but the outer por- 
tions are causally linked to the inner regions through angular 
momentum transport, and it is these outlying regions that set 
the long-term accretion disk time scales. 

The physical agent responsible for angular momentum 
transport and turbulent dissipation of energy is now believed 
to be the magneto-rotational instability (MRI) or Balbus- 
Hawley instability (Balbus & Hawley 1998) which involves 
the shearing amplification of a weak seed magnetic field. The 
genera] process of angular momentum transport and energy 
dissipation in accretion disks is still commonly referred to as 
“viscosity”. Balbus & Papaloizou (1999) show a foundation 
for the dynamical “a" view of accretion disks (Shakura & 
Sunyaev 1973) within the MRI. 

Extensive study of accretion disks in interacting binaries 
has guided thought on accretion disks in other types of sys- 
tems, such as around BHs in AGNs. For instance, the best 
constraint we have of the magnitude of the a parameter comes 
from a detailed study of the rate of decay of dwarf nova out- 
bursts (Smak 1984), and the values obtained for ionized gas 
a ~ 0.1 — 0.2 have been carried over to accretion disks in 
other systems. 

One aspect of accretion disks which is quite different for 
isolated BHs is the lack of a constraint on the outer radius. 
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The sudden introduction of an annulus of gas near a central 
mass leads to a process of smearing of the matter into a broad- 
ened ring. The ring spreads both to smaller and larger radii in 
response to outward transport of angular momentum. As mass 
drains from the inner edge, an ever smaller amount of mate- 
rial moves to larger and larger radii. Pringle (1974) studied 
“external accretion disks” in which matter initially in a torus 
spreads both to smaller and larger radii. The class of solu- 
tions he considered of relevance to this work is that for disks 
with constant total angular momentum and decreasing mass. 
Pringle (1991) generalized his earlier results. 

By writing the equations for mass continuity and angular 
momentum transport in cylindrical coordinates, assuming Ke- 
plerian rotation ft 2 K = GA/bii?’~ 3 and a thin accretion disk, 
and integrating over the vertical thickness of the accretion 
disk, one arrives at an equation for the evolution of the sur- 
face density E = 2 ph, where p is the density and h the disk 
semithickness (actually pressure scale height). 


3E __ 3 d_ 
dt r dr 



y>i / 2 



(0 


The kinematic viscosity coefficient 

2aP 


' 3fW (2) 

where P is the pressure and a is the Shakura-Sunyaev 
parametrization of the angular momentum transport and heat- 
ing (Shakura & Sunyaev 1973). Debate persisted over many 
years as whether the viscosity should scale with the total pres- 
sure or gas pressure. A scaling with the total pressure leads 
to a strong dynamical instability in regions where P l3 .d > 
P gas — the Lightman-Eardley instability (Lightman & Eard- 
ley 1974) — which manifests itself as a limit cycle operat- 
ing in the inner, radiation pressure dominated portions of the 
disk (e.g., Cannizzo 1997). However, the fact that X-ray bina- 
ries with very high rates of accretion, such as Sco X- 1 , often 
have perfectly stable light curves leads one to conclude that 
the Lightman-Eardley instability is probably the artifact of an 
incorrect assumption for thin disks, namely v oc P to tai (Done, 
Gierliriski, & Kubota 2007). With the advent of the MRI, de- 
tailed MHD calculations show that the saturation limit of the 
instability depends on the magnetic energy density reaching 
a limit that may depend on surface density rather than pres- 
sure (e.g., Hirose, Krolik, & Blaes 2009). Hirose et al. find 
not only that radiation dominated disks are thermally stable, 
but that fluctuations in magnetic field strength causally pre- 
cede those in pressure (both radiation and gas). Hirose et al. 
suggest adopting E as the relevant independent variable. 

The MRI couples linearly to the Keplerian shear in the disk 
and dissipates energy through turbulent transport. In formu- 
lating the steady state equations for the radial profile of accre- 
tion disk structure, one begins by equating a heating function 


with a cooling function 
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assuming an optically thick vertical structure with energy 
transported by radiative diffusion. The vertically integrated 
stress provides the total radiated flux (assuming no radial ad- 
vection of energy) 
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Finally, the condition of vertical hydrostatic equilibrium in the 
disk gives 

^ = n 2 K h 2 . ( 6 ) 

p 

2.1 . Blandford-Znajek Coupling to the Inner Disk 

After a period of equilibration and expansion, the exter- 
nal disk may span many decades in radius — from ~ 10 7 
to ~ 10 11 cm. The physical state of the disk is expected to 
vary enormously as one progresses from small to large radii. 
In terms of mediating the flow of mass down onto the cen- 
tral engine, the similarity solutions just discussed rely on the 
physical state of the outer disk. The physical state of the in- 
ner disk determines the strength of the Blandford-Znajek (BZ) 
process (Blandford & Znajek 1977) that extracts rotational 
energy from the central BH. Ghosh & Abramowicz (1997) 
examined the dependence of the strength of the BZ mecha- 
nism on physical conditions in the inner disk. The energy 
density of the magnetic field entrained within the disk sets 
the BZ luminosity, and insofar as the saturation limit is set by 
equipartition with the pressure, the maximum pressure p m&x 
is of relevance. 

Livio, Ogilvie, & Pringle (1999) investigated the role of 
the BZ mechanism in powering outflows and found that if the 
magnetic field threading the inner disk does not differ signifi- 
cantly in strength from that threading the black hole, the elec- 
tromagnetic output from the inner disk regions should dom- 
inate over that from the hole. In simple terms, the disk is 
more efficient at holding and accelerating the magnetic field 
because its magnetic diffusivity is much less than that of the 
BH. 

However, more recent works using sophisticated GRMHD 
codes to follow the detailed magnetohydrodynamic launch- 
ing of flows from the very inner disk find them to be highly 
efficient in terms of tapping the available supply of rest mass- 
energy in the disk. (McKinney & Narayan 2007ab; Krolik, 
Hawley, & Hirose 2007). McKinney & Narayan (2007a) note 
that previous simple estimates of the BH and disk failed to 
take into account the conversion of EM power into thermal 
and materia] power in the corona and disk wind (Ghosh & 
Abramowicz 1997, Livio, Ogilvie, & Pringle 1999) and so 
may have seriously overestimated the power from the disk, 
in which case the BH power may dominate. Krolik, Haw- 
ley, & Hirose (2007, see their Table 1) find high efficiencies 
r i, earn for the jet power, measured either in terms of mass, 
radiative energy, or EM energy, normalized to the rest mass- 
energy accretion in the disk Me 2 . Their radiative efficiency 
values (which they denote t?nt) increase with a/M, ranging 
from 0.057 for a/M — 0 up to 0.264 for a/M = 0.99. The 
efficiencies may increase with time due to BH spin-up. These 
workers find the jet to be Poynting flux dominated near the 
BH spin axis, giving strength to the idea of a lightweight jet 
that responds quickly to changes in M in the inner disk. In 
this work we adopt £i, ean . = 0.1 as a fiducial value relating 
the jet energy production as seen in X-rays to the accretion 
rest mass-energy within the main body of the accretion disk. 

Initially the accretion rate can approach ~ 1M© s _i , 
i.e., super-Eddington by a factor ~ 10 14 . One might ques- 
tion whether an accretion disk description is fully applicable. 
Studies of the inner disk at these very early times find that the 
temperatures and densities are so large that the cooling is by 
neutrino emission (Narayan, Piran, & Kumar 2001, Kohri & 
Mineshige 2002, Kohri, Narayan, & Piran 2005, Chen & Be- 



4 


CANNIZZO & GEHRELS 


loborodov 2007). Chen & Beloborodov (2007) show that for 
a BH with high spin a ~ 1 a neutrino cooling dominated disk 
transitions to an advective disk at ~ 10~ 3 — 10 2 Mq s -1 . 
Kohri & Mineshige (2002, see their Fig. 1) detail the differ- 
ent physical regimes of these disks. The cooling function 

Q~ — Grad + Qadv + Qv ’ ( 7 ) 

where Q~ is the neutrino cooling. 

2.2. Eddington Luminosity and Accretion Rate 

For the GRB depicted schematically in Fig. 1,2 = 0.54 
and dh ~ 3.1 Gpc. The inferred, beaming-corrected X- 
ray luminosity corresponding to the earliest time we seek to 
address, namely the onset of segment II, would be L\\ — 
4tt d 2 /beam Pxrt. ii - 3 x 10 43 erg s~\ where we used 
Pxrt, ii — 10- 11 erg cm" 2 s~ : (Gmpe et al. 2007; see 
their Fig. 4) and /beam = 3x 10~ 3 . The isotropic-equivalent 
X-ray plateau luminosity for this GRB, ~ 10 46 erg s _1 , is 
in about the middle (logarithmically) of the range 10 44 erg 
s~ 1 < 47rd|_FxRT, ii < 10 48 erg s~ J shown in Dainotti et al. 
(2008; see their Figs. 1 and 2). Adopting fbeam = 0.1, the 
geometry-corrected luminosity translates into a rate of accre- 
tion in the inner disk L \\/ ~ 5.3 x 10~ 3 M© yr -1 . 
In the conventional view, for rates of accretion larger than the 
Eddington rate 

Lb — ~ 1.25 x 10 39 (Mbh/ 10M©) ergs -1 

K'es 

( 8 ) 

radiation pressure impedes the accretion process (see below, 
however, for the possibility of a slim disk). The Edding- 
ton accretion rate depends on radius, through the local effi- 
ciency with which gravitational energy is radiated e acc (r) = 
GMbh/ [ re 2 ). The locally defined Eddington accretion rate 
Me{t) is related to the Eddington luminosity via 

L e = G.c c(r)M E (r)c 2 . (9) 

Thus e-s.ee ( r ) <x r -1 and M e ( r) ex r, varying from 
1.3 x 10 ~'Mq yr _1 in the inner disk 6GMun/c 2 to 1.5 x 
10~ 3 M© yr -1 at 10 11 cm. Therefore in the region of interest 
to us in this study, namely the outermost part of the disk, the 
rate would only be locally super-Eddington by a factor ~ 3 at 
this early time. 

The preceding discussion was based on the standard as- 
sumption of solar composition material, which is almost cer- 
tainly not the case here. The fail-back matter comprising the 
disk would be expected to be dominated by r — processed ma- 
terial - Ni, Si, and Fe group elements. Calculations of the 
Eddington limit for such material indicate a potential suppres- 
sion of the standard Eddington value by up to about three or- 
ders of magnitude due to increased opacity from lines (e.g.. 
Fryer, Colgate, & Pinto 1999). However, this applies mainly 
to much lower temperatures than those discussed in the next 
section. 


2.3. Steady State Physics 
2.3.1. Standard Thin Disk 

Cannizzo & Reiff (1992) present general Shakura-Sunyaev 
scalings for disks with arbitrary opacity law. For optically 
thick disks for which v cx P ga6 , k = K es = 0.2(1 + A') — 0.2 
cm 2 g -1 (no hydrogen), and mean molecular weight p = 16 
(appropriate for the expected for mixture of alpha elements 


- predominantly 56 N, 28 Si, 16 0, and 12 C - Fryer, Colgate, & 
Pinto 1 999), the surface density 

S = 2.84 x 10 7 g cm -2 m mi i a ~l° (10) 

where r n = r/10 11 cm, itibh, i = Mbu/10Mq, a_i = 
a/0.1, rh— 2 = M/10~ 2 Mq yr -1 . (In the region of interest 
for this study rra _ 2 — 0.01 — 1.) The midplane temperature 

T = 7.83 x 10 6 K r~ 9/w a 1 1/5 ml’/. (11) 


The disk temperature is less than the virial temperature 

T viri ai = \ £ — = 4.26 x 10 8 K r// m BH , 1 • (12) 

6 it r 


The only assumptions needed up to this point are Q~ = 
Q~ = Q. To obtain the disk scale height and density, we 
now specify P — P gas . This yields a pressure scale height to 
disk radius ratio 


h 

r 


0.055 r\[ 2Q m. 


7/20 -1/10 . 1/5 


'BII, 1 a -l 
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-2 > 


(13) 


a density 


p = 2.57 x 10 gem r 


-3 ^-33/20 11/20 —7/10 .2/5 


11 
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BII, 1 a -l 


m. 


■2 1 
(14) 


and self-gravity to central object gravity 
il _ = 1.62 x IQ" 3 rT° m-// 2 ? a/ 7 / 10 ^ 


9c 

The radiation pressure to gas pressure ratio 


m_ 2 
(15) 


Prad 


21/20 7/20 1/10 . 4/5 


- gas 


— 91 r n ^bpi, 1 a -i 


m_ 


(16) 


exceeds unity, therefore the pressure scale height, which was 
predicated upon P = P ga s. should be increased by a factor 
192 « 10 to h/r ~ 0.5 because of hydrostatic equilib- 
rium. Note that because we adopt 12 oc P gas , the fact that 
Prad > Pgas is not problematic in itself as long as h/r « 1. 
However, the fact that h/r ~ 0.5 begins to be problematic. 

As a last consistency check, the ratio of the flux of radially 
advected energy (M /27rr 2 )c 2 to vertically transported energy 


Qadv _ 3 K es 1l MS 

Qrad 16 tt acp r 2 T 3 

1 ao ia — 3 1/10 — 7 / 10 — 1/5 * 2/0 / 10 \ 

= 1.02 x 10 r n m BH l a_ l m_ 2 . (18) 


In summary, the standard model is marginally self-consistent 
for representative values in the outer disk rn ~ 777.-2 — 1- 
Note that for ?nBii. 1 = rn = 1, we have e ac0 (r) ~ 
1.5 x 10 -5 and Mb{t) — 1.5 x 10 _3 Mq yr _! . A value 
777—a — 0.5 (determined from Lu) at rn = 1 would be 
locally super-Eddington only by a factor ~ 3, but by larger 
factors at smaller radii since Afjs(r) oc r. As noted ear- 
lier, the true super-Eddington factors could be greater for the 
1 — processed gas, although for T ~ 10 7 K the opacity is basi- 
cally electron scattering, which is actually suppressed due to 
the fact that the hydrogen fraction X ~ 0, and K es oc (1 4 - X). 
At smaller radii the standard model would become quite prob- 
lematic due to Prad/Pgas >> 1 leading to h/r > 1, and 

M/Me{t ) » 1. Although a detailed description of the disk 
at smaller radii is not of primary relevance to this work, we 
briefly look at two formalisms of potential interest. 
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2.3.2. Super-Eddington Slim Disk 

For disks with accretion rates of order Eddington, the siz- 
able radial gradients demand the introduction of a radial ad- 
vective energy transport term (Abramowicz et al. 1988, Chen 
et al. 1995; Kato, Fukue, & Mineshige 1998). The cooling 
function is now written 

= + O 9 ) 

where the advective transport 


and s is the entropy per particle. This can be written as 

QZl V =^>r- p C, (2D 

where £ is a logarithmic derivative involving s that is of order 
unity (Chen et al. 1995), or 

M 0 

Qadv ~ 2 1Tr 2 c «’ ^ 2 ) 

where c 2 = (3/2 )GM/r. The anticipated accretion rate ~ 

3M E (r) at rn = 1 places the accretion disk model on the 
upper stable branch of slim disk solutions (Abramowicz et al. 
1988, see their Fig. 4) where radiation pressure dominates 
(Kato, Fukue, & Mineshige 1998, see their section 10.1.1). 
Solving the disk parameters for P — P rac i gives 

E = 9^0^ (23) 

— 3.87 x 10 4 g cm' 2 m B aZ\ m_ 2 , (24) 


3 

4n arot 


1.23 x 10 G K r x 1 :/S tobh, i 1//8 rhi^ 4 , 


0.86G. 


, n -7 _3 — 3/2 — 1/2 —1 . 

p — 2.23 x 10 g cm r n m BH 3 a_ x m— 2 , 
§21 = 4.05 x 10 s r;* /s m ^ 


(Shakura & Sunyaev 1973), the maintenance of a disk border- 
ing on being Eddington at all radii would enforce h/r — 1. 
The disk would be non-steady state, 

32 7t c 


and strong winds would be generated at all radii. Solving for 
the disk variables with P — P ra ,i and h/r — 1 yields 

S - — ^ — (33) 

3aUA' K es r 


■ 6.94 x 10 4 gem 2 r\ x m 


BII, 1 u -l> 


T =(l~ y - 


1.16 x 10 6 K r u 3,/S mBii, i 1,/S 


p = 3.47 x 10 7 g cm 3 r u 1/2 m B ^ 2 a_}, 


= 2.20 x 10 3 r xx '° m B jj , a. 


Since there is no dependence on accretion rate, this state 
would persist, in theory, until M dropped below M E at the 
outer edge. An annular region with constant M(r) would 
then begin to propagate to smaller radii. This appears at 
first glance to raise the intriguing possibility of applying 
this long period of constant accretion onto the central en- 
gine to explain the plateau phase of the light curve shown 
in Fig. 1 (segment II). If each annulus in the disk emits 
at roughly its local Eddington limit, an integration over the 
disk shows that the total accretion luminosity can exceed L e 
by a factor ~ ]n(A/[r outer ]/A'f£[nnner]) (Shakura & Sun- 
yaev 1973; Begelman, King, & Pringle 2006) which would 
have a value ~ 10 in the current application. However, the 
inferred GRB segment II luminosity L n — 3 x 10 43 erg 
s -1 (/beam/3 x 10~ 3 ) greatly exceeds even this enhanced Ed- 
dington value ~ 10 x Z, E (10M©) ~ 10 4D erg s~\ unless 
the beaming factor were orders of magnitude smaller than the 
value 3 x 10~ 3 we take in this study. Such a tightly collimated 
(low would not be expected from an Eddington outflow. 


Q a j v _ ^ 1 M 

Qriul 127tcr 3 Af E (r) 


1.12 r,/ rn —O' 


Thus the starting assumptions P ra< ] > P sas and Q a(h , > Q rad 
are validated. 

2.3.3. Marginally Eddington Disk with P = P rad 

In view’ of the super-Eddington accretion rates, another pos- 
sibility worth investigating is a disk in which strong mass out- 
flow' due to radiation pressure driven winds occurs as the local 
accretion rate begins to exceed M E . Since the aspect ratio 

< 31 , 

’■ WeH 


2.4. Time Dependent Behavior of External Disks 

Pringle (1974, 1991) presented exact self-similar solutions 
for the time dependent behavior of disks w'ith a free outer 
boundary for the viscosity parameterization 


Of interest are the class of solutions for which the disk total 
angular momentum is constant. 

2.4. 1 . Standard Model 

For the simplest case - gas pressure and electron scattering 
opacity — one may combine eqns. (2)-(4) to find 

v ~ CrTr- 3 . (41) 
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where 


C 


v 4 / 3 


3 ( k l 


4/3 


(- 


\ 


1/3 


V vrrip ) V 1 2acGM m] J 


(42) 


This is in the form considered by Pringle (1991) for m = 2/3 
and n = 1. The solution is 


nr-t) 

So 




(43) 


expanding exp(— x) as 1 — x and combining the t dependen- 
cies gives a steeper decay, Lj oc t~ 7 / 3 . The dividing point 
^divide = fo/9 between L d oc f~ 4 / 3 and La cx t~ 7 ! 3 can be 
evaluated by using the closure relation between C^, A), an d 
to. yielding 

^divide — = 5 /0 S mgj^ J , (53) 

where flo — Qk( t o)- 


where 


2.4.3. Irradiation-dominated Disk 


f{u) = (28)~ 3/2 u~ 3/h (l - v 7/5 ) 3/2 (44) 

(Cannizzo, Lee, & Goodman 1990), which can be verified by 
substitution into eqn. (1). The dimensional scalings £ 0 , r 0 , 
and to satisfy 

fo-roC-^o 273 . (45) 

The solution given in eqn. (43) extends out to u = 1, so that 
the spreading of the outer edge varies as (r/ro) — {t/to) 3 ' 8 . 
The right hand side of eqn. (43) multiplies the £ solution 
from eqn. (10). The total disk angular momentum is constant, 
and the disk mass decreases as 

4 / t \ “ 3 / 16 

M d (t) = (28)- 3 / 2 yr 0 2 £ 0 (^J , (46) 


At very late times corresponding to segment IV, which 
is not always seen, one needs also to consider irradiation- 
dominated disks. For disks in which the locally defined irra- 
diation temperature T irr = [Lacc/^tror 2 )] 1 / 4 , where L acc = 
£ acc Me 2 , exceeds the disk midplane temperature, the viscous 
heating of the disk will be subordinate to irradiational heating. 
One now finds 4 

v = C lrr r 4/3 E 1/3 , (54) 


where 


Co 


a 


,4/3 


k B V /3 ( 4e acc c 2 
H ni p J \ 27 oG 2 


1/3 


(55) 


This is in the form considered by Pringle (1991) for m — 1/3 
and n = 4/3. The solution is 


and therefore the available accretion luminosity L,i(t) oc 
Ma oc f- 19 / 16 . 


2.4.2. Super-Eddington Slim Disk 
In this model 

v = C &< wr 1/2 t (47) 

where 

Cadv = . (48) 

Form = 0 and n = 1/2 the solution from Pringle (1991; see 
his eqn. [3.2.8]) is 


nr,t) 

Eo 




- 4/3 

(i -ku b y. 


(49) 


where u — (r/ro) 1/,2 (#/io) _1,/3 . k — m(4m + 4 — 2 n) -1 
(5m + 4 — 2n)~ 1 oc m = 0, a — 1/m — + oc, and 6 — 3. 
Euler’s formula 


lim (l+— — e x (50) 

M-OC fij 

enables us to recast £ as 

IV 1 ]. 
to J | 

(51) 

Integration of 2trnirE(r, f) for the disk mass gives 




(52) 


For t « t o /9, the exp term will go to zero. Hence one has 
A/,./ oc f~ 1/3 and therefore La oc in agreement with 

Kumar et al. (2008a) for s — 0. For late times t >> to/9, 



where 

f{u) = (18)- 3 u- 2 (l-t i 2 ) 3 . (57) 

The disk mass 

M d {t) = { 18)- 3 |r 2 E 0 ^) V3 , (58) 

so the disk luminosity Ld oc t~ 4 ! 3 . 

3. GENERAL SCENARIO 

We argue that the entire X-ray sequence spanning > 10 6 
s can be understood in terms of the external accretion disk 
passing through a sequence of different physical states. We 
propose that the jet becomes active almost immediately after 
creation of the BH and the beginning of formation of the fall- 
back disk. The subsequent decrease in X-ray luminosity is 
then due to the decreasing mass supply in the inner disk. The 
current paradigm is that a jet of material launched from the 
core of the progenitor decelerates as it encounters circumstel- 
lar material (Sari et al. 1999; Frail et al. 2001). 

In the following discussion of the four Zhang et al. regimes, 
we adopt GRB 060729 as a specific example due to its very 
long coverage both in terms of flux decay and time (both span 
~ 7 decades — see Grupe et al. 2007 for a detailed study). 
The light curve (Grupe et al. 2007, see their Fig. 4; Will- 
ingale et al. 2007, see their Fig. 10) shows a 0.3 - 10 keV 
flux initially at ~ 3 x 10~ 7 erg cm“ 2 s _I at ~ 100s, de- 
creasing to ~ 2 x 10 -11 at ~ 10 3 s, then down slightly to 
~ 10~ n by ~ 5 x 10 4 s, and finally to ~ 3 x lO - " 14 at 

4 Menou, Pema, & Hemquist (200 1 , see their Appendix A2) neglect the M 
(and hence E) dependence of u, and therefore effectively derive a similarity 
solution for a disk with constant irradiation. 
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~ 5 x 10 6 s. The redshift 2 = 0.54 gives, for the stan- 
dard cosmology, a luminosity distance di =3.1 Gpc, so that 
4-?r d\ ~ 1.15 x 10 5 ' cm 2 . To take into account cosmologi- 
cal time dilation the long term light curve depicted schemat- 
ically in Fig. 1 should also be contracted along the x - axis 
by a factor 1 + z = 1.54. Adopting a geometric beam- 
ing factor /beam ~ 3 x 10~ 3 , a beam efficiency ebeam = 
i'beam /(Me 2 ) ~ 0.1, and neglecting the bolometric correc- 
tion, the plateau flux Fxrt, ii — 10~ u erg cm" 2 s'" 1 trans- 
lates into an isotropic luminosity 47rd 2 Fxrt, ii = 1.2 x 10 lfi 
erg s -1 , a beaming corrected luminosity Lu — 3.6 x 10 43 
erg s _1 , and a rest mass-energy flux in the disk Me 2 — 
L n/fbeam = 3.6 x 10 44 erg s'" 1 , which would not be di- 
rectly observable presumably due to the advective nature of 
the inner disk. This value is significantly in excess of the Ed- 
dington value Le e; 1.3 x 10 39 (Mbh/ 10M®) erg s _1 . It is 
also important to acknowledge that the very long duration of 
the plateau forGRB 060729, ~ 50, 000s, is not typical. Dain- 
otti et al. (2008) find a range between about 10 2 s and 10 5 s in 
plateau duration, corrected to the source frame. 

Taking /beam = 3 x 10~~ 3 and eteam = 0.1 for entire 
evolution shown in Fig. 1 (which is probably unrealistic), 
the amount of matter accreted during each segment is 1.5 x 
1O“ 4 M 0 for the prompt emission out to 100s, 5.5 x 10 -5 M© 
for segment I, 8.5 x 10“ 6 M® for segment II, 2.4 x 10 " 5 M q 
for segment III, and 1.7 x 10~ fi Af® for segment IV. It is im- 
portant to note that these numbers represent the bookkeeping 
associated with the powering of the jet, which we are assum- 
ing produces the X-rays. Given the likelihood of significant 
mass depletion within the disk, e.g., M(r) oc r s , with s in 
the range 0 to 1 (Kumar et al. 2008a), the actual total mass 
available from fail-back could be orders of magnitude larger, 
~ 0.1 — 1M®. Given the long term persistence of the external 
disk, the logarithmic rate of decay for each of the four seg- 
ments identified by Zhang et al. (2006) must be relatable to 
different physical regimes the disk passes through during its 
long period of mass depletion: 

3. 1 . Segment I: Fall-back Debris 

a\ ce 3 for 10 2 < t < 10 3 s: During this early time the 
inner disk is being dynamically assembled out of fall-back 
material, and so the assumptions entering into the similar- 
ity solution for the evolution of the external disk are not yet 
valid. Kumar, Narayan, & Johnson (2008a) considered the 
evolution of debris from the SN event that disrupts the pro- 
genitor and creates the BH and subsequent GRB, and show 
analytically that the arrival of material onto the BH from 
the collapse of the outer half of the progenitor core leads to 
M ~ r 3 . In their model the rapid infall of matter arises be- 
cause of the steep density profile of the outer core. During 
this early period of extremely high accretion rate > 10 10 £e 
initially, the inner disk may contain regions where the cooling 
is neutrino-dominated and then advection-dominated. During 
the early, chaotic phase of fail-back and accretion, the high 
rate of accretion significantly depletes the amount of material 
surrounding the BH. 

It has been suggested that early in the evolution of the disk 
self-gravity may play a role (e.g., Perna, Armitage, & Zhang 
2006). Therefore it may be possible to account for this early 
evolution purely within the confines of an accretion disk for- 
malism, through considering self-gravitating disks (J. Pringle, 
2008, private communication), although it is not clear a priori 
whether such a steep decay law can be obtained. If there is 


still substantia] material at small radii left over from the im- 
mediate post-GRB accretion (MacFadyen & Woosley 1 999), 
it may be possible to model this decay time as a disk with self- 
gravity driving the viscosity. Detailed numerical calculations 
of self-gravitating disks support the contention that the self- 
gravity torque can be efficient enough to maintain the disk 
at the edge of instability to self-gravity, and that the relevant 
physical length scales are small, of order the disk thickness h, 
in which case a local viscosity formalism can be utilized (e.g., 
Lodato & Rice 2004). 


3.2. Segment II: Transient Plateau 
an ~ 0.5 for 10 3 < t, < 10 4 s: The isotropic luminos- 
ity observed during the plateau 10 46 erg s" 1 translates into 
a mass accretion rate of ~ 5 x 10~ 3 Mq yr -1 after apply- 
ing the corrections given earlier. We argue that the plateau 
stage occurs when a small amount of material forms a torus 
at the radius of the progenitor, ~ 10 11 cm. The similarity 
solutions discussed earlier apply to accretion disks that have 
had time to adjust into a quasi-equilibrium profile which is 
close to steady-state. Therefore in practice there will be a pe- 
riod of adjustment during which time the small amount of gas 
(> 3 x 10~ 5 Mo) that once was part of the envelope of the 
progenitor becomes stretched out due Keplerian shear into an 
accretion disk. The time dependent calculations of Cannizzo 
et al. (1990) indicate a flat transient period in the light curve 
in response to an annulus of material being suddenly intro- 
duced near a central mass. An apparent problem with the 
light curves shown in Cannizzo et al. (1990, see their Fig. 
3) is that the transient plateaus shown there have slightly in- 
creasing fluxes, whereas the Swift XRT plateaus are flat or 
decaying. The initial conditions for their time dependent cal- 
culations were tori with physically ad hoc Gaussian radial 
profiles, E(r) oc exp[-(r - r c ) 2 / (2w) 2 }, with w/r c — 0.5. 
Therefore time was required after the calculations started for 
matter to spread significantly to smaller radii, approach the 
standard disk profile £(r) oc r" 3 ' /4 , and raise the central ac- 
cretion rate onto the BH. For the present problem, however, 
it is more likely that one begins with significant amounts of 
gas at smaller radii. Therefore the accretion rate onto the BH 
would reach a higher level more quickly, and the light curves 
would be flat or decaying. Obviously, time-dependent ac- 
cretion disk calculations for the present problem are desired, 
were a physically motivated initial £(r, v) profile available 
from SN fail-back calculations. 

Since our scenario begins to come into play at the beginning 
of segment II, the accreted masses quoted earlier for segments 
II through IV must add up to a minimum for our initial disk 
mass, i.e„ about 3.5 x lO -5 !/®. The relation between £o. r o< 
and t 0 given in eqn. (45), in combination with the numerical 
normalization from Cannizzo et al. (1990, see their Fig. 3), 
can be used to form a scaling for this interval 


to 


5 x 10 4 


c r 7 / 3 n 4 / i! \r 
br 0. Ii a -l M d 


1/3 

'BH, 


(59) 


where r 0 , a = r o /10 n cm, and A/,?. - 4.5 = A/u.sk/10"' 4 5 

J/. : . We argue that this plateau period corresponds to Zhang 
et al’s phase II. 


3.3. Segment III 
3.3.1. Standard Disk 

am — 1.2 for 10 4 < t < 10 5 s: This segment corresponds 
to the standard model for the outer disk for which the solu- 
tion given earlier applies. Gas left over from the progenitor 
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envelope at ~ 10 11 cm has become smeared out and evolved 
into the radial profile of an accretion disk, and now the clas- 
sical similarity solution for the decay begins to apply. The 
disk mass varies as tr 3j/16 so the luminosity, which scales 
with d,M,i /dt, varies as f~ 19 /' 1G ~ t~ 12 . For other opacity 
laws k — K 0 p a T b , Cannizzo, Lee, & Goodman (1990) found 
d log L/d log t = - (38 + 18a - 46)/(32 + 1 7 a - 26). This re- 
duces to the electron scattering limit (—19/10) fora = 6 = 0. 
In the optically thick limit, the decay rate is relatively insen- 
sitive to the opacity law. For instance, a Kramer’s law a — 1, 
b — —3.5, gives a decay rate d log L/d\og t = -5/4. 

3.3.2. Super-Eddington Slim Disk/ADAF-like state 

For GRBs with high rates of accretion in the post fall- 
back disk, conditions may be so extreme that the standard 
model is not relevant. Therefore the super-Eddington slim 
disk must be considered. In section 2.4.2 we present the 
time dependent solution for such disks, and show that for 
early times d\ogL/d\ogt — —4/3, which would probably 
be indistinguishable in the XRT data from the standard disk 
value d,\ogL/d\ogt = —19/16. Kumar, Narayan, & John- 
son (2008a; see their sections 3.1 and 3.2) take a different 
approach in calculating the time dependent evolution of the 
transient disk. Rather than solving the E(r, t) evolution equa- 
tion they prescribe a functional form M(r) oc r s , and then 
solve for the total disk angular momentum, mass, and outer 
radius. The motivation for this M (r) form is that for “ADAF- 
like” conditions one expects strong mass outflow driven by a 
positive Bernoulli constant (Narayan & Yi 1994, 1995). Thus 
,s = 0 would represent no mass loss, and s — 1 would corre- 
spond to a flow in which 90% of the mass is lost from the disk 
for each inward decade in radius. This technique bypasses the 
need for a detailed treatment of the E(r, t) evolution equa- 
tion and the kinematic viscosity v. It also ignores any po- 
tential transient phase through which the disk might have to 
evolve, and therefore assumes that the disk goes immediately 
into its equilibrium state. They derive a luminosity decay law 
Lj(t) oc f“ 4 ( s+1 )/ 3 , which reduces to t~ 4//3 for s = 0. This 
agrees with our law for super-Eddington slim disks at early 
times. 

3.4. Segment IV 
3.4.1 . Irradiation Dominated Disk 

qiv — 2 for t > 10 6 s: Menou, Perna, & Hernquist (2001) 
consider the evolution of supernova fail-back disks and find 
that irradiation becomes important at late times t > 10 6 s. 
Since the radial decrease in irradiation temperature is flatter 
than that for disk temperature (r -1 / 2 versus r~ 3 / 4 ), as the 
disk mass decreases and the surface densities and tempera- 
tures in the outer disk continue to drop, at some point the con- 
dition T\ tt > Tviscous will be satisfied, and the standard model 
similarity solution will be supplanted by the irradiation model 
solution. Moreover, a strong wind will begin to deplete sig- 
nificantly the disk mass, at a rate that exceeds the local mass 
flow (Begelman, McKee, & Shields 1983, Begelman & Mc- 
Kee 1983, Shields et al. 1986). Therefore the decay rate will 
steepen from the standard model value a ~ 1.2. Begelman, 
McKee, & Shields (1983, see their eqn. [4.2]), derive a char- 
acteristic mass loss rate due to a Compton heated wind 




where £ is a collection of factors of order unity, and L cr ~ 
0.03Te- Taking a luminosity L — eMc 2 gives 

M c h - 2 x lO -3 M© yr -1 (e_im_ 2 ) 2/3 n 4 h, n ( 6i ) 

which will serve as a strong source of mass depletion in the 
outer disk. 

The time-dependent decay solution for the wind-driven 
phase presented in section 2.4.3 is not complete because it 
does not include the effect of mass decrease on E(r, t). Also, 
the decay law, etjv = 4/3, is the same as the early-time so- 
lution for the super-Eddington slim disk, therefore could not 
account for the steepening decay for segment IV. The gener- 
alized Kumar et al. (2008a) decay solution in the presence of 
a disk wind, L<i{t) oc t~ 4 ( 1+a )/ 3 t provides more insight into 
the potential effect of a wind. The fact that when a segment 
IV is observed it has a steeper slope than segment III sug- 
gests that if there are ADAF-like conditions in the disk early 
on, mass loss is not significant since the observed decay rate 
agrees with theoretical one if a ~ 0. Since the ADAF-like 
state is supposed to be driven by a positive Bernoulli con- 
stant, as the disk mass drains and physical conditions within 
the disk become less extreme, one would not expect for mass 
loss to become stronger if the positive Bernoulli constant were 
the only force driving mass loss. However, in the picture of 
Menou. Perna, & Hernquist (2001), in which irradiation be- 
comes important at later times as a strong agent in driving a 
wind, a power law m(r) oc r s to characterize the variation of 
mass flow within the disk in the presence of strong evapora- 
tion may become relevant. The steeper decay law in segment 
IV would then indicate a significant non-zero s value; for ex- 
ample a nominal value s ~ 0.5 would be required to obtain 
Ld(t) oc 

3.4.2. Super-Eddington Slim Disk, late time solution 

In section 2.4.2 we found that for the time dependent super- 
Eddington slim disk, the disk mass fori > O.l(afio) -1 varies 

as Ma(t) oc t~ 4 / 3 , and therefore Lj(t) oc t ~ 7 / 3 . This is close 
to the aiv — 2 from Zhang et al. (2006) and suggests an- 
other possibility for the late time decay. In section 2.3.1 we 
presented the standard scalings for usual accretion disk with 
v cc P gas and k = K es . A look at the various self-consistency 
tests showed that most of them were at least marginally vali- 
dated for conditions expected in the outer disk during segment 
II. However, the range of isotropic-equivalent luminosities for 
segment II (Dainotti et al. 2008) can range up to ~ 2 orders 
of magnitude greater than that adopted in section 2.3.1 , which 
would place an uncomfortable strain on the standard model. 
For GRBs with outer disks at such extreme conditions that the 
standard model does not apply, the super-Eddington slim disk 
maybe a better description. The range in L\\ values shown 
in Dainotti et al. (2008) appear to show an intrinsic scatter, 
therefore there may be GRBs for which the standard model 
does apply at late times (low L\\), and others for which the 
slim disk model is required (high L\\). One prediction of this 
idea would be that intrinsically brighter bursts would be more 
likely to exhibit a segment IV, with its steeper decay. How- 
ever, for the fainter bursts one would not be able to follow the 
late-time decay to as low a flux level, and therefore it might 
be difficult to make a definitive observational statement about 
the presence of a segment IV. 

3.5. Mass supply: Viscous accretion vs. Fail-back 


Mch — 4.5 x 10 ' A/© yr ' 


1/3 
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Kumar et al. (2008b) present and discuss the criteria for de- 
termining whether viscous accretion or fail-back will be the 
primary supply mechanism feeding mass into the central en- 
gine. Both processes will influence the light curve, but the 
one with the longer time scale is likely to have the dominant 
effect. After the SN has occurred, the fall-back time for an el- 
ement of gas in the progenitor envelope to fall from radius r to 
the BH is roughly the free-fall time £fb ~ 2(r 3 /GMBn) 1 ^ 2 - 
(This ignores any residual pressure support during the col- 
lapse, and other effects such as a wind from the accretion 
disk and a shock originating near the BH during its forma- 
tion.) Inverting this expression gives the radius in the star 
expected to be accreting at a given time £ after the SN oc- 
curs, rio — . , where £2 = i/100 s. The spe- 

cific angular momentum j*,i8 = 3. 6m B jj 1 r 1 ' 0 /n(r), where 
s is in units of 10 18 cm 2 s^ 1 , and fn(r) is the local angu- 
lar velocity of a gas at radius r in the progenitor, in terms of 
the Keplerian value fl#-. After the SN explosion, gas which 
ends up being trapped in the potential well of the BH will 
fall inward to the point where it hits its angular momentum 
barrier, determined by fa(r). This viscous accretion time 
scale £ acc ^ (r-//i) 2 (an^[r circ ])“ 1 « 2o“ I OA'' 1 (r circ ) 
(since h/r — 1), where r c i rc is the circularization radius 
~ j1/(GM\m). This line of reasoning leads to eqn. [4] of 
Kumar et al. (2008b), 

facc/ffb ^ 10 aZlfn 3 (r). (62) 

Kumar et al. (2008b) note that this relation implies the accre- 
tion time will be shorter than the fall-back time, and therefore 
subordinate as a mass supply mechanism, for fn(r) < 0.4. 
They find that the steep decay segment (I) would be possible 
only if the disk can adjust rapidly to a decrease in A/ft,, and 
use this consideration to place an upper limit on the rotation 
speed in the core (on the assumption that segment I is due to 
fail-back of the progenitor core). Their approach is strength- 
ened by our finding that all plausible disk states give decay 
rates with power laws < 2, too shallow to explain segment I. 

At early times before the disk has had a significant oppor- 
tunity to spread, the fundamental physical considerations pre- 
sented by Kumar et al. (2008b) seem well-motivated. At later 
times, however, the increasing viscous time scale associated 
with the outer disk edge, and the maintenance of a causal 
coupling to smaller radii, via angular momentum transport, 
begin to introduce a significant correction factor that must 
be taken into account in eqn. (62). It is important to note 
that implicit in the derivation of eqn. (62) is that the radius 
appearing in the fail-back time and the viscous time are the 
same. However, the spread of the outer disk edge introduces a 
strong nonlocality: instead of taking £ acc oc a~ l O.K~ l \r c irc ) 
one should more properly consider o/~ 1 QK~ A {'^disk, outer)- 
Since £ acc varies as r 3/2 , this introduces a correction fac- 
tor (r c usk, outer/t’circ) 3 '' 2 into eqn. (62). Taking a fiducial 
spreading rate (r/r 0 ) = (i/£o) 3 / 8 as in the standard solu- 
tion, and adopting typical numbers for times separating the 
different segments, £j/n 10 3 s, tn/m — 5 x 10 4 s, and 
fin/iv - 3x 10 6 s, this implies an increase over the locally 
defined (i.e., at r cjrc ) accretion time scale by a factor ~ 8 at 
£ii/m and ~ 90 at £iji/iv> respectively (stemming from in- 
creases of -ruisk, outer by a factor 4 from ty jj to £ij/m, and 
by a factor 20 from tyw to fm/iv)- The faster rate of expan- 
sion calculated by Kumar et al. (2008a), (■ r/r 0 ) = (£/£o) 2/3 , 
would give larger correction factors. Also, these factors are 


actually lower limits because they only take into account the 
expansion of the disk relative to £jyn- The disk expansion en- 
hances the effect of the accretion supply with respect to the 
mass supply by prolonging £ acc , and also facilitates an under- 
standing of why the steepest decay (segment I) must occur 
at the earliest time possible, before the significant outward 
movement in r^-sk, outer has enforced £ acc >> £ft,. 

Kumar et al. (2008b) present examples of 4 Swift GRBs 
for which ani — 4 — 6, much steeper than in the standard 
Zhang etal. (2006) picture am ~ 1.2. For these rare bursts 
the Kumar et al. fail-back scenario may well provide a better 
description than our viscous accretion disk scenario. 

4. DISCUSSION 

The overarching question behind this study is the follow- 
ing: Is the very long term (£ > 10 e s) decay seen in the X-ray 
flux from GRBs due to the deceleration of a relativistic jet, 
or to a decreasing mass supply onto the central engine? We 
argue in favor of the latter. We have presented a general sce- 
nario for understanding the long term evolution, in particular 
the four power law decay segments found through Swift XRT 
observations. The basic idea relies on the fact that the ever- 
expanding outer edge maintains a causal connection through 
the intervening radii all the way to the inner disk, given the 
coupling of the viscosity mechanism to the local shear in the 
disk. In more basic terms, the mass of the accretion disk de- 
creases as a power law in time, resulting in a decrease in mass 
accretion rate at the inner edge. 

The overall evolution is as follows: 

(I) The initial steep decay phase arises from the dynam- 
ical fall-back of ~ 0.1 M© of material into the vicinity of 
the BH (Kumar, Narayan, & Johnson 2008ab). Kumar et al. 
2008b show by comparing the viscous and fall-back limes in 
the early stages of fail-back that £ acc << £fb- The disk there- 
fore responds quickly to the steeply decreasing, externally im- 
posed mass supply from the progenitor core. This scenario is 
only possible at early times before expansion of the outer ac- 
cretion disk causes £ acc > £fb. An alternative but perhaps 
less likely possibility for the steep decay involves completely 
different physics (Pe’er et al. 5 2006). 

(II) The initial establishment of the disk leads to a brief 
transient phase during which calculations indicate a flat light 
curve. This corresponds to the viscous smearing and spread- 
ing of gas left over from the progenitor envelope into an ac- 
cretion disk (Cannizzo et al. 1990), which we take to lie at 
roughly the progenitor radius ~ 10 11 cm. The dynamical time 
scale there QjT ~ 10 3 s is toward the lower range of observed 
plateau durations (Dainotti et al. 2008, se their Figs. 1 and 2), 
and our scaling for the transient period following the arrival 
of material to establish the accretion disk is in line with the 
duration of the plateau phase seen in some of the XRT light 
curves. 

The plateau phase has been ascribed by previous workers 
to continued activity of the central engine, or to the accretion 
disk having low value of the Shakura-Sunyaev a parameter 
(e.g., Kumar, Narayan, & Johnson 2008a; although in their 
final analysis they downgrade this possibility). Given the ap- 
parent robustness of the non-linear saturation limit of the vis- 
cosity parameter in the MRI process (Beckwith, Hawley, & 

1 Their explanation, which involves a hot cocoon with modest Lorentz 
factor and a nearly thermal spectrum, may be refuted by Swift observations 
showing that the GRB prompt emission joins smoothly to segment I — an un- 
likely coincidence given that the hot cocoon represents an independent phys- 
ical component. 
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Krolik 2008), however, we also view the small a explanation 
as unlikely. Kumar, Narayan, & Johnson (2008a) view as a 
more plausible explanation continued fall-back onto the cen- 
tral BH. Kumar, Narayan, & Johnson (2008b; see their Fig. 2) 
present conditions for the core/envelope density structure of 
the progenitor that would be required to obtain the plateau (as 
well as the steeper segments - 1 and III). In the previous sec- 
tion we showed how the spreading of the outer disk enhances 
the effective viscous accretion time, relative to the fall-back 
time, and ensures that viscous accretion will become increas- 
ingly dominant as the outer disk expands significantly beyond 
what had been the radius of the progenitor. 

(III) After the transient period, the outer disk becomes de- 
scribable using the standard formalism (t~ 1,2 decay) or an 
advective formalism (f -1 - 3 decay). The continued high rate 
of accretion substantially depletes the mass of the disk. The 
inner disk is still not only substantially super-Eddington (by 
a factor ~ 10° — 10 7 ) but also potentially degenerate due to 
the very high densities (e.g., Narayan, Piran, & Kumar 2001, 
Kohri & Mineshige 2002, Kohri, Narayan, & Piran 2005, 
Chen & Beloborodov 2007, Ohsuga & Mineshige 2007). 

(IV) At some point the outer disk may become not only ra- 
diation dominated, but also subject to a substantial Compton 
driven wind, leading to an increase in the decay rate. If the 
rate of mass flow within the body of the disk can be charac- 
terized M(r) oc r s as in Kumar et al. (2008a), then their 
decay solution La(t) oc f- 4 ( 1 +®)/ 3 would provide a good fit 
for 8 ~ 0.5 regardless of whether the disk was standard or 
advective. Alternatively, our late time super-Eddington slim- 
disk solution L,i(t) oc t ~ 7 ' 3 may apply at this time. 

Depending on details such as the initial mass of the disk, 
states (II) and (IV) may or may not become manifest in the 
long term light curve. 

The current paradigm for the long term GRB light curves 
attributes changes seen in the decay properties to variations 
(i.e., a decrease) in the bulk Lorentz flow of a jet of material 
expelled from the central engine as it collides with and inter- 
acts with circumstellar material. We argue that a strong jet is 
established almost immediately after the GRB prompt emis- 
sion, and the subsequent decay properties of the light curve 
are set by a dwindling mass supply due to the long term de- 
crease in E at all radii. The decay law obtained by Kumar et 
al. (2008a) in which Lx oc where the mass loss 

rate from the disk varies as r s , allows for a range of decay be- 
tween t~ 4 ' 3 and t~ 8 / 3 . The steeper decay, as an explanation 
for segment III, seems to be excluded by the data (e.g.. Ghis- 
ellini et al. 2009, see their Fig. 16). If the slim disk/ADAF 
picture proves to be a valid description of the outer disk, the 
observed segment III decay rate would appear to argue against 
a strongly varying mass loss rate with radius. Also, the poten- 
tial bimodality in decay rate with time found in section 2.4.2 
for the slim disk state, t~ 4/3 for early times and for late 
times, might play a role in the III/IV transition. The dividing 

■> 3/2 __ i / 2 _ 1 

point lies at ~ 10''s r 0 | n r?? BH ' ja_j, so for this to work the 
disk would have to have expanded to > 10 12 cm to account for 
the long time scale ~ 10 J — I0 r> s of the III/IV transition. Al- 
ternatively, the onset of a strong wind at late times, driven by 
irradiation, may force s to increase from 0 to ~ 0.5, thereby 
steepening the decay in segment IV. 

In addition to our alternative explanation for the plateau 
seen in some of the XRT light curves, there is also an is- 
sue raised by the high-latitude emission model (Kumar & 


Panaitescu 2000). In this model, after the mechanism respon- 
sible for producing X-ray emission ceases, an observer along 
the central jet axis sees continued emission as photons from 
increasingly off-axis lines of sight continue to arrive. For uni- 
form surface brightness, an intrinsic spectrum produces 
a power law decay in time t~ a , where a — /? + 2. This pre- 
diction was tested by O’Brien et al. (2006, see their Fig. 4) 
for the initial decay segment (Zhang et al’s I) using a large 
sample of Swift GRBs. Instead of a ~ (3 4 - 2, they found 
a scatterplot in a versus (3. Thus the standard model also 
seems problematic for decay segment I as regards the tem- 
poral/spectral evolution, although more recent work on the 
high-latitude “curvature effect” has been promising (Zhang et 
al. 2007,2009). 

5. IMPLICATIONS 

Eqn. (1) of Sari, Piran, & Halpem (1999) has proven of 
great use to observers insofar as it relates observable quan- 
tities, namely jet break time, redshift, and GRB fluence, to 
a beaming angle (assuming a density of the circumstellar 
medium). The outline of the theory we present is not yet 
sufficiently developed to build a complex formalism. Simple 
statements can still be made, however. Our eqn. (59) gives an 
approximate scaling for the time of the segment II/III transi- 
tion. If the luminosity of the plateau scales roughly with the 
mass accumulated during fail-back that supplies the transient 
disk associated with the plateau (segment II) in our model, 
then eqn. (59) predicts an inverse scaling between L x ,n and 

<n, in, namely Lx, 11 oc f n ni . (This ignores the r 0 scaling, 
however.) An inverse correlation is observed in the data, but 
appears less steep (L x ,u oc t jj°f n ) than our prediction (Dain- 
otti, Cardone, & Capozziello 2008; see their Figs. 1 and 2). 

The Kumar et al. (2008a) decay law with mass loss from 
the disk, La(t) oc is consistent with our early 

time decay law for the super-Eddington slim disk, and with 
the observed value am ~ 1.2, for s = 0. It is also close 
to the decay law for standard disks, L,i(t) oc t~ 19 / 16 . This 
appears to indicate that for very high M disks where the stan- 
dard model is not applicable during segments II and III, if an 
ADAF-like state is correct, then mass loss must be minimal 
(i.e., s ~ 0, or M(r ) ~ constant). The steepening then as- 
sociated with segment IV could either be due to the onset of 
significant evaporation, e.g., s ~ 0.5 would be sufficient, or 
the late time super-Eddington slim disk law Ld{t) oc 

If our model has some merit, then a revision in current 
thinking would be needed. First, there may not be a stan- 
dard GRB energy reservoir as in Frail et al. (2001). Also, 
for the jet to respond quickly to whatever conditions prevail 
at the central engine (i.e., accretion rate), it may need to be a 
light-weight jet with a low degree of baryon loading. It may 
be a Poynting-flux dominated jet as in the model of Lyutikov 
& Blandford (2003). The deceleration experienced by such a 
jet would tend to be less, for a given energy flux. The pre- 
vious estimates of the circumstellar density field n(r) would 
need to be revised. Conversely, the jet may be baryon-loaded 
during the steep decay segment (segment I) which appears to 
be too steep to be covered within an accretion disk formalism. 
The spectral evolution observed during this phase may be dif- 
ficult to account for with a purely EM jet (O’Brien et al. 2006; 
Pe’er et al. 2006). There may also be some baryon loading 
during the rest of the decay, but the density of the circumstel- 
lar medium would need to be less than previously thought, so 
that deceleration would be minimal. 



In this work we have implicitly assumed that the X-ray light 
curve represents the bulk of the radiated jet energy. There- 
fore our model runs into the same difficulty as the standard 
model in accounting for why the breaks seen at other wave- 
lengths do not follow the X-ray. The physical emission mech- 
anisms must differ across different wavebands, which could 
also potentially mitigate one of the criticisms of the standard 
model. One idea currently being discussed is that, within a 
given GRB there can co-exist jets with different opening an- 
gles and different physical conditions, each of which can con- 
tribute preferentially to emission in different wavebands. If 
the opening angles are considerably different, this could lead 
to different break times (e.g., Racusin et al„ 2008). Our sim- 
ple model is not yet developed enough to make any predic- 
tions on the relative rates of decay at different wavelengths. 

There are other important issues to consider, such as the 
formation and stability of an accretion disk in the middle of a 
supernova. Enough material would have to be evacuated from 
the vicinity of the disk near the BH at early times to permit the 
disk to survive the expected shocks and radiation of the explo- 
sion. Also, if the BH acquires a kick velocity r^ick during the 
initial explosion, this velocity would have to be small enough 
so that the accretion disk was not left behind. (A simple con- 
sideration of the binding energy shows that matter in the disk 

with n or bit = 1155 km s" 1 x rj‘ 1 1/2 > u kick , where 
t>kick — 100 - 300 km s~\ should remain bound to the BH.) 
Finally, the jet issuing from the BH would have to be power- 
ful enough to punch through the veil of material comprising 
the ejecta, and to maintain an open channel for > 1 0 6 s. If the 
jet is primarily electromagnetic as opposed to baryon-loaded, 
the open EM field lines may help in maintaining a clear path 
for the jet beam. 

Recently various workers have proposed a simple way for 
obtaining the flat decay associated with segment II. Implicit 
in our work and that of others such as Kumar et al. (2008ab) 
is that the start of the prompt emission is the relevant start 
time for all the other emission - X-ray, optical, and radio. It 
is possible that the physical mechanism responsible for pro- 
ducing the X-rays actually commences ~ 10 4 — 10 5 s before 
the prompt emission, so that the flat appearance of segment 
II is merely due to having an incorrect start time (Yamazaki 
2009; Liang et al. 2009). Yamazaki (2009) plots four ex- 
amples of GRBs for which a shifting of the start time for the 
X-ray decay produces a good fit for the plateau, after having 
optimized the decay for segment III. One can however find 


other examples of GRBs for which there is a gradual decay 
in segment II that is not well-fit by simply shifting the start 
time, therefore this explanation may not work for all GRBs. 
Also, it remains unclear as to how the X-ray emission could 
begin significantly before the GRB, which presumably marks 
the moment of explosion of the progenitor and formation of 
the BH. 


6. CONCLUSION 

We present a new paradigm for understanding the long-term 
behavior of GRBs, as inferred primarily through Swift XRT 
observations. We argue that the four power law segments into 
which the long term XRT light curves can be divided can be 
explained by the time dependent evolution of the external ac- 
cretion disk. The initial, steep period of decay happens before 
the external accretion disk has been established and therefore 
is beyond our scope (although a self-gravitating disk remains 
a possibility). Simple estimates of the accretion and fail-back 
time scales face and ifb indicate that fall-back will be dom- 
inant at the earliest times. Therefore the fact that segment I 
has the steepest slope, too steep to be accounted for by the 
disk models we study, strengthens the picture advocated in 
Kumar et al. (2008b) in which the direct fall-back of the pro- 
genitor core leads to segment I. After the initial evolution, 
expansion of the outer disk edge enforces t acc > so that 
accretion is dominant. We ascribe the subsequent flattening 
in the GRB light curves to a period of viscous re-adjustment 
in the accretion disk formed out of the gas elements from the 
progenitor envelope that remain bound to the BH and end up 
at their specific angular momentum radii after the hypernova 
explosion. The following steepening corresponds to the sim- 
ilarity solution for the decay of a standard disk or advective 
disk at the largest radii, and the further slight steepening at 
later times could be due to the additional mass depletion from 
a strong wind, or in some cases to the late time evolution of a 
super-Eddington slim disk. 
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